LS I +61 303 as a potential neutrino source on 
the light of MAGIC results 



Diego F. Torres 1 ' 2 & Francis Halzen 3 
February 5, 2008 

Institucio Catalana de Recerca i Estudis Avancats (ICREA) 

Institut de Ciencies de l'Espai (IEEC-CSIC), Facultat de Ciencies, Universitat Autdnoma de 
■arcelona, Torre C5 Parell, 2a planta, 08193 Barcelona, Spain. Email: dtorres@ieec.uab.es 
Department of Physics, University of Wisconsin, Madison, WI 53706 

Abstract 

Very high energy y-rays have recently been detected from the microquasar LS I 
+61 303 using the MAGIC telescope. A phenomenological study on the concomi- 
tant neutrinos that would be radiated if the y-ray emission is hadronic in origin is 
herein presented. Neutrino oscillations are considered, and the expected number 
of events in a km-scale detector such as ICECUBE is computed under different 
assumptions including orbital periodicity and modulation, as well as different pre- 
cision in the modeling of the detector. We argue that the upper limits already 
imposed on the neutrino emission of LS I +61 303 using AMANDA-II and the 
forthcoming measurements by ICECUBE may significantly constrain -in an inde- 
pendent and unbiased way- the y-ray to neutrino flux ratio, and thus the possibility 
of a hadronic origin of the y-rays. The viability of hadronic models based on 
wind-jet interactions in the LS + 61 303 system after MAGIC measurements is 
discussed. 

1 Introduction 

The possibility of an hadronic origin of high energy radiation from microquasars 
and gamma-ray binaries have been recently extensively discussed (e.g., see Romero 
et al. 2001, 2003, Bosch-Ramon et al. 2005 and references therein). Earlier works 
have already dealt with neutrino and high energy emission from galactic sources, 
and especially those presenting jets. Among them, Levinson and Blanford (1996) 
where among the first to suggest a possible correlation between microquasars and 
gamma-ray sources. Later, Levinson and Waxman (2001) showed that if the en- 
ergy content of the jets was dominated by an electron-proton plasma, a several hour 
outburst of 1-100 TeV neutrinos produced by photomeson interactions should pre- 
cede the radio flares that are associated with major ejections. Photopion production 
in jets was also considered by Distefano et al. (2002), by assuming the parame- 
ters of their model as inferred from radio observations. A recent review on the 
high energy aspects of astrophysical jets, together with more complete references 
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to previous works, can be found in the work by Levinson (2006). Partially, recent 
work has been motivated by the discovery of the presence of relativistic hadrons in 
microquasar jets like those of SS 433, as inferred from iron X-ray line observations 
(e.g. Migliari et al. 2002). In addition, some sources (as LS 5039 and LS I +61 
303) present quite inconspicuous behavior at X-rays, with low variability and flux 
levels. Protons, generally being subject to less efficient energy loss mechanisms 
than electrons, could be accelerated to high energies generating only low luminos- 
ity counterparts at other frequencies. Microquasars have then been proposed to be 
an additional source of cosmic rays (e.g., Heinz & Sunayev 2002). It is however 
true that the proof of protons being behind the highest energy photons detected 
so far has yet to be given at a level that not even reasonable doubt is tenable. At 
the moment, neither the analysis of EGRET detections correlated with supernova 
remnants (see Torres et al. 2003 for a review) nor the recent HESS results (e.g., 
Aharonian et al. 2006) are conclusive in this respect. This is a central problem 
of cosmic-ray physics. However, key multiwavelength observations (for instance, 
for RX J0852. 0-4622, see Aharonian et al. 2006b) are not far from elucidating the 
origin of the gamma radiation, at least for some cases. It is clear that the detection 
of concomitant neutrino emission from a hadronic source of gamma-rays would be 
a definite proof. 

Here, under the model independent assumption that the y-ray spectrum mea- 
sured by the Major Atmospheric Gamma Imaging Cherenkov telescope (MAGIC, 
Albert et al. 2006) results from inelastic pp collisions, a lower limit to the neu- 
trino yield expected from LS I +61 303 is computed. Disregarding the nature 
of the compact object in LS I +61 303, recently classified as a gamma-ray pulsar 
wind system by Dhawan et al. (2006), a hadronic interpretation of the gamma 
emission would result in a prediction of a neutrino yield. In any case, if there are 
jets, their matter content is unknown. We note, however, that modeling the high 
energy radiation with pp collisions is not discarding leptonic contributions (see 
e.g., Bosch-Ramon et al. 2006), since those might dominate at lower gamma-ray 
energies and at other locations in the orbit of the binary system, where electron 
losses are diminished. The appeal of hadronic emission of TeV radiation in close 
X-ray binary systems can be understood if one considers the losses and acceler- 
ation process for electrons, under the assumptions of a given magnetic field. We 
discuss these issues in more detail below. In this paper, we discuss how this pp 
assumption can be tested using neutrino telescopes under different scenarios for 
variability, and whether earlier proposed hadronic microquasar models can still 
comply with current observational constraints. The key approach we propose here 
is that an experiment such as ICECUBE could test the level of enhancement of 
neutrino emission as compared with that detected at TeV y-ray energies. By com- 
paring with theoretical modeling of photon absorption, ICECUBE can ultimately 
test or disprove a hadronic scenario for LS I +61 303. 

2 LS I +61 303: basic facts 

LS I +61 303 shares with LS 5039 the quality of being the only two known mi- 
croquasars that are spatially coincident with sources above 100 MeV listed in the 
Third Energetic Gamma-Ray Experiment (EGRET) catalog (Hartman et al. 1999), 
and the only two detected so far at higher y-ray energies (Aharonian et al. 2005, 
Albert et al. 2006). These sources both show low X-ray emission and variability. 
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Optical spectroscopic observations of LS I +61 303, which is located at a distance 
of 2.0+0.2 kpc (Frail & Hjellming 1991), revealed a rapidly rotating BO main se- 
quence star with a stable shell and radial velocities compatible with binary motion 
in a 26.5 day orbital period (Hutchings & Crampton 1981). The most accurate 
value of the orbital period, P orb =26.4960±0.0028 d, comes from the analysis of 
more than 20 years of periodic radio outbursts (Gregory et al. 2002). The maxi- 
mum of the radio outbursts varies between phase 0.45 and 0.95 (assuming T =1D 
2.443.366,775) following a modulation of 4.6 years. An X-ray outburst starting 
around phase 0.4 and lasting up to phase 0.6 has also been detected (Goldoni & 
Mereghetti 1995, Taylor et al. 1996, Harrison et al. 2000). The periastron takes 
place at phase 0.23 and the eccentricity is 0.72+0.15. Extended jet-like and ap- 
parently precessing radio emitting structures at angular extensions of 0.01-0.05 
arcsec have been reported by Massi et al. (2001, 2004); this discovery has sup- 
ported the microquasar interpretation of LS I +61 303. The uncertainty as to what 
kind of compact object, a black hole or a neutron star, is part of the system (e.g., 
Casares et al. 2005), seems to have been lifted while this paper was under review 
by Dhawan et al. (2006). These authors have presented preliminary results from a 
July 2006 VLBI campaign in which rapid changes are seen in the orientation of a 
cometary tail at periastron. This tail is consistent with it being the result of a pulsar 
wind (in which case, the gamma-emission should be the result of the interaction of 
this wind with the companion star outflow). No large features or higher velocities 
were noted on any of the observing days, i.e., there was no notice of the high 0.6c 
velocity outflow as reported by Massi et al. 2005, which implies at the least its 
non-permanent nature. The changes within 3 hours were found to be insignificant, 
so the velocity can not be much over 0.05 c. In fact, tail velocities around 7500 
km s _1 were measured at periastron. If confirmed, these results would classify LS 
1+61 303 as a pulsar wind system (a gamma-ray binary) instead of a microquasar 
(powered by mass accretion). Debate is alive. 

3 MAGIC detection of LS I +61 303 

The spectrum derived from MAGIC data between 200 GeV and 4 TeV at orbital 
phases between 0.4 and 0.7 is fitted by a power law function: 

F y = (2.7 ± 0.4 ± 0.8) x 10- 12 (£/TeVr 2 - 6±0 - 2±0 - 2 cnrV'TeV-', (1) 

with the errors quoted being statistical and systematic, respectively (Albert et al. 
2006). MAGIC measurements showed that the very high energy y-ray emission 
from LS I +61 303 is variable. The maximum flux corresponded to about 16% of 
that of the Crab Nebula, and was detected around phase 0.6 with 8.7<x of signifi- 
cance. In addition, since the system was observed in different orbital periods and 
the detections occurred at similar orbital phases, there is a hint (although not yet 
a published proof) of a periodic nature of the VHE y-ray emission. A definitive 
proof or disproof of periodicity of the y-ray emission in LS I +61 303 awaits the 
publication of further MAGIC observations. 
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4 Timescales and models 



When analyzing the possibility of a leptonic origin of the gamma-ray emission (in 
such a case the neutrino production would be quenched), the location of the emit- 
ter within the source and the dichotomy between a fast jet or a slower cometary 
tail wind are crucial. To understand why, one can compare the relevant timescales 
and impose constraints upon the magnetic field. The maximum energy of acceler- 
ated electrons is generally computed by requesting that the acceleration timescale 
equals that of the losses (a.k.a. cooling), where the latter is provided by the 
electron's radiation via synchrotron and Compton mechanisms. The acceleration 
timescale is (e.g., Malkov & Dury 2001, Aharonian et al. 2005) 

W = n ^ * 0. 1 1 £rev B a l rj s, (2) 

where r L = E/eB is the Larmor radius, E TeV = Ej 1 TeV, and B a = B/ 1 G is 
the strength of the ambient magnetic field, rj gives account of the efficiency of 
acceleration: in the case of extreme accelerators (maximum possible acceleration 
rate allowed by classical electrodynamics) r\ — » 1, whereas for shock acceleration 
in the Bohm diffusion regime r\ ss 10(v/c)~ 2 (see e.g. Malkov & Dury 2001). For 
the Compton scattering we consider both, the Thompson and the Klein-Nishina 
regimes. In the former (say, for E e <<c 0.1 TeV) the cooling time is inversely 
proportional to the electron energy, 

fxhom * 0.03£ Tc ' v s. (3) 

In the latter, the characteristic Compton cooling timescale can be approximated by 

t KN * 34wo'£?ev s > ( 4 ) 

where wo = w r /500 erg cm 4 . The Compton scattering of TeV electrons against 
starlight occurs in the Klein-Nishina regime. By equating the acceleration and 
Compton timescales, r acc = r K N> we obtain 

E e , max - 82657.9[B G (v/c)V] 33 TeV - 2[fi G (v/0.2c) 2 w '] 3 3 , (5) 

where the first equality is obtained just replacing the previous formulae and do- 
ing the algebra, and the second equality provides a useful scaling especially for 
microquasar jets, where velocities of 0.2c are tenable (e.g., as in LS 5039). The 
companion star of the compact object in LS I +61 303 is a BO V with a tempera- 
ture about 22500 K. Its optical luminosity is taken to be L+ » 2 x 10 38 erg/s, and 
periastron and apastron are reached at about 2.5 and 14.5 stellar radii (e.g., Romero 
et al. 2005 and references therein). Assuming that R+ =s \QR a , as appropriate for 
a Bo V/Be star, periastron (apastron) occurs at 1.7 x 10 12 cm (1.0 x 10 13 cm). The 
energy density of the starlight 

varies then between 200 and 5 erg cm 4 . Formally, then, even for slow winds, a 
magnetic field B > 1 G allows for the maximum energy of accelerated electrons to 
approach/exceed 10 TeV. For such fields however, synchrotron emission also plays 
a role. The synchrotron timescale is 

t sync * 400S G 2 £ Tc ' v s, (7) 
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that is, for instance, at Bq ~ 10, the cooling is dominated by synchrotron radiation 
rather than Compton processes. Together with the equality f acc = f sync , gives 

£e.™x * 19S G l/2 (v/c)TeV. (8) 

These equations imply that electrons could be accelerated to energies exceeding 10 
TeV in a low-magnetic field environment and/or low stellar density w . In the case 
of a jet system, low magnetic fields can be attained far from the inner parts of the 
jets, out of the core of the binary system. A periodic feature in the gamma-radiation 
(which is evidence of y-ray production inside the binary system) would then likely 
favor hadrons, which cool less efficiently, as the primary particles. Their maximum 
energy would be determined by the condition r L < R^, which gives 

E p < 3 X 10 15 (tf Jct /10 8 cm)(S/10 5 G) eV, (9) 

i.e. and inner-jet magnetic field of about 10 5 G can produce PeV primaries. These 
protons are not sufficiently energetic to interact with starlight to trigger photome- 
son processes, and since the X-ray emission from LS I +61 303 is not indicative of 
the existence of an important accretion disk, the photomeson process against the 
X-ray photons is expected to be sub-dominant as well (the X-ray luminosity is only 
about 10 33 erg s -1 ). It is also possible to see that pion production would inevitably 
lead to neutrino emission (otherwise stated, that pions decay before interacting in 
the media), and also that the production of neutrinos from the subsequent muon 
decay also proceed with high probability as long as the magnetic field does not 
exceed B ~ 10 6 G (see the computations from Aharonian et al. 2005 for the case 
of LS 5039: since LS I +61 303 have similar X-ray fluxes, they can be adapted vis 
a vis to this scenario). 

Following the recent results form Dhawan et al. (2006), LS I +61 303 would 
be similar to the system PSR 1259-63/SS2883. In the latter, a 48 ms radio pulsar 
is found in an eccentric orbit around a 10th magnitude main sequence star with an 
strong wind, a Be star with T cSf ~ 27000 K. LS I +61 303 would have, however, 
a much shorter period. The PSR 1259-63 system has been observed by CANGA- 
ROO (Kawachi et al. 2004) and HESS (Aharonian et al. 2005) and modeled by 
both, a mainly leptonic (e.g., Kirk et al. 1999) and a mainly hadronic contribution 
(Kawachi et al. 2004). Interestingly,the observed lightcurve of the PSR 1259-63 
system seems to be qualitatively similar to the prediction made by the model of 
Kawachi et al. (2004) where hadronic interactions and neutral pion production in 
the mis-aligned stellar disk plays a dominant role in the gamma-ray production 
mechanism. Neronov et al. (2006) have also presented a hadronic interpretation of 
the PSR 1259-63 HESS data, where the TeV emission is the result of pp collisions 
between protons accelerated in the pulsar wind and nuclei resident in the wind of 
the stellar companion. Albeit a full hadronic model of LS I +61 303, in the case 
it actually is a pulsar wind system, is still missing and is beyond the scope of the 
present work, one certain sub-product of a hadronic gamma-ray production would 
be the concomitant emission of neutrinos. In what follows, we analyze in a model 
independent way -starting from the MAGIC data- up to what level a detector like 
ICECUBE could detect such neutrinosQ 

1 NEMO/KM3NET, having a planned effective area similar in size to ICECUBE -see, e.g., Distefano 
et al. 2006 for a general discussion- would also be a key player in the possible detection of this system. 
However, due to location in the southern hemisphere, ICECUBE is better suited for an observational study 
on LS I +6 1 303 and we focus the discussion on it. These alternative instruments would of course be key for 
the detection of southern hemisphere sources, such as LS 5039. 
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5 Neutrino yields 



5.1 Neutrino fluxes 

Inelastic pp collisions lead to charged (2/3) and neutral (1/3) pions. Thus, if 
hadronic interaction are involved in the production of the y-ray flux, a compa- 
rable emission of and v e is to be expected. For a detailed discussion of the 
computations of pp collision y-rays see, for example, the appendices in the works 
by Torres (2004) and Domingo-Santaman'a & Torres (2005). The neutrino flux at 
the production site can be obtained from the y-ray flux measured at Earth, although 
with the caveat that the latter can be affected by absorption. Thus, the so-predicted 
neutrino flux is to be considered as a lower limit to the one actually produced. In 
the GeV and TeV regime, if y-rays are significantly absorbed, the y-ray spectrum 
may be steeper than the neutrino one, what further secures the lower limit qual- 
ity of the estimation of the neutrino flux when starting from that in y-rays. The 
v-fiux is F' v = k'Fy where k' = (Q' v /Q y ) stands for the ratio of -at the source- 
emissivities of the particular kind of neutrinos i = p, e, r or their antineutrinos, and 
y-rays. This ratio depends on the value of the photon spectral index, which (before 
absorption of y-rays proceeds) is to be shared by the neutrino spectrum. The most 
recent derivation of k' has been presented by Cavassinni et al. (2006), using the 
numerical package PYTHIA. Details of the prescription we use below and the way 
in which it was derived can be found there. Galactic distances are much larger 
than the neutrino oscillation length for energies around TeV. Then, flavor oscilla- 
tion probabilities have to be considered in computing the neutrino flux at Earth. 
After propagation in vacuum, the original neutrino flux will be modified according 
to Fj = Y,j=e,!i,T wnere F°. is the expected neutrino flux in the absence of 

oscillations, and Pij the probabilities of inter-conversion among neutrino flavors 
(Cavassinni et al. 2006). Assuming that the CP violating phase is negligible (e.g., 
Constantini and Vissani 2005), oscillations almost completely isotropize the sig- 
nal, and generate a non-negligible v T flux at Earth, even when it is negligible at the 
source. 

Taking into account the measured (average) MAGIC spectrum of y-rays, the 
predicted neutrino fluxes at the source out of this spectrum, following the previ- 
ously commented derivation, are 
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whereas the fluxes at Earth after taking into account v-oscillations are 
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Errors in the y-ray spectrum (about 30% in the normalization and 10% in the slope, 
including statistics and systematics in quadrature) would propagate directly into 
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the neutrino flux. It is to be stressed that the latter neutrino fluxes are lower limits 
only. Since absorption of y-rays in the stellar photon fields of the binary system 
is unavoidable, in the case of a hadronic production of y-rays, the neutrino flux 
ought to be larger than the y-ray one (e.g., Aharonian et al. 2005, Christiansen et 
al. 2006). In the case of LS 5039, for instance, Aharonian et al. (2005) studied 
the potential of such source as a neutrino emitter after the measurement by HESS, 
and proposed that the neutrino flux expected above 1 TeV could be up to a hundred 
times larger than the y-ray flux detected in that energy band. In such an increased 
range, ANTARES, an experiment with similar effective area for neutrino detection 
than AMANDA (see below), would be able to test the putative neutrino emission. 
In the case of LS I +61 303, and in the framework of the hadronic microquasar 
model developed by Romero et al. (2003) and Torres et al. (2005), Christiansen et 
al. (2006) reported that the opacity-corrected y-ray flux above 350 GeV dropped 
from 7.1 x 10~ n cm~ 2 s" 1 to 1.4 x 10 -12 cm~ 2 s _1 , at the periastron of the system, 
i.e., a reduction by a factor of 50 in y-ray flux, that would not affect the associated 
neutrino production (we discuss further this model below). Neutrino telescopes 
can constraint this very factor of enhancement, the neutrino to photon ratio. 



5.2 Testing emission scenarios with neutrinos 

Indeed, we argue that the combination of the y-ray flux measured by MAGIC from 
LS I +61 303 with the existing upper limits from neutrino experiments already 
restricts the enhancement factor, and thus provides constraints for a hadronic origin 
of the y-ray radiation from LS I +61 303. We first note that the neutrino flux of 
Eq. i ll It is consistent with AMANDA-II upper limits on this source (Ackermann 
et al. 2005). The latter has been imposed at 90% CL as 2.4 x 10~ 8 cm" 2 s~\ 
integrated for energies above 10 GeV and assuming an spectral slope of -2. The 
quoted limit correspond to the 2000-2002 data sample. This upper limit is not 
very restrictive, although it is only a factor of 25 larger (disregarding the possible 
difference that the slope used to obtain it may introduce, which does not affect 
the concept) than the flux given in Eq. (2). Ackerman has recently presented a 
new limit from 4 years of AMANDA II data, at a level of 6 x 10~ 9 cirr 2 s' 1 . The 
combined F V/l and F v lower limit flux yielded by LS I +61 303 at Earth as deduced 
in Eq. Jilt , when integrated above 10 GeV, is only about a factor of 6 smaller than 
the improved AMANDA-II upper limit. Then, we can compute how much time the 
next generation of neutrino experiments would need to constrain further a possible 
hadronic emission; i.e., how well can neutrino experiments constrain the neutrino 
to photon ratio0 To do this, we briefly summarize how to approximately compute 
the signal and background events for an ICECUBE-like detector. 



5.2.1 Neutrino detection 

Neutrino telescopes search for up-going muons produced deep in the Earth, and 
are mainly sensitive to the incoming flux of v M and v^. ICECUBE will consist of 

2 If the neutrino spectrum is harder than the photon one (as assumed for instance in the case of LS 5039 
by Aharonian et al. 2005), the neutrino flux would be even larger than just the scaling up of the y-ray flux, 
so that to impose a constraint over the neutrino to y-ray ratio, the most conservative choice is that the change 
in spectrum is not significant, which is an assumption we follow. For instance, if the increased neutrino flux 
is due to absorption, the neutrino flux would be similar to the one we use but with a standard E~ 2 spectrum 
(see, e.g., Halzen & Hooper 2005). This will increase the number of events. 
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Table 1 : Computed neutrino yields for two years of observations in ICECUBE assum- 
ing the neutrino to photon ratio is equal to 1 . 



Orbital scenario 


Signal Events 


Background Events 


Strict Periodicity 


0.6 


5.0 


Modulation 


1.3 


10.0 



4800 photomultipliers, arranged on 80 strings placed at depths between 1400 and 
2400 m under the South Pole ice (e.g., Halzen 2006). The strings will be located 
in a regular space grid covering a surface area of 1 km 2 . Each string will have 
60 optical modules (OM) spaced 17 m apart. The number of OMs which have 
seen at least one photon (from Cerenkov radiation produced by the muon which 
resulted from the interaction of the incoming v in the earth and ice crust) is called 
the channel multiplicity, iV ch . The multiplicity threshold is set to N cb = 10, which 
corresponds to an energy threshold of 200 GeV, the same threshold of the MAGIC 
y-ray observations. The angular resolution of ICECUBE will be ~ 0.7°. A first 
estimation of the event rate of the atmospheric y-background that will be detected 
in the search bin is given by (e.g., Anchordoqui et al. 2003, Romero & Torres 
2003) 

— b = A eff J dE v P V ^(E V ) Afl , (12) 

where A eff is the effective area of the detector, Afl « 1.5 x 10~ 4 sr is the angular size 
of the search bin, and d® B /dE v < 0.2(£ r /GeV)- 3 21 GeV" 1 cirr 2 s" 1 sr 1 is the 
Vf, + atmospheric v-flux (Volkova 1980, Lipari 1993). Here, P r ^(E v ) denotes 
the probability that a v of energy E v on a trajectory through the detector, produces 
a muon. For E v ~ 1 - 10 3 GeV, this probability is « 3.3 x 10~ 13 (£ v /GeV) 2 " 2 , 
whereas for E v > 1 TeV, P V ^(E V ) a 1.3 X 10~ 6 (£ r /TeV) 08 (Gaisser et al. 1995). 
On the other hand, the y-signal is similarly obtained as 



dN 
dt 



= A eff J dE v (F Vfl + F fll ) P y ^(E v 



) , (13) 



where (F v + F 9/i ) is the incoming v^-flux. In the previous integrals we use both 
expressions for Py^(E v ) according to the energy, and integrate from 200 GeV up 
to 10 TeV. The approximate validity of these expressions for the probability can 
be indirectly verified by the comparison of the neutrino background estimates and 
measurements in experiments such as ICECUBE. 



5.2.2 Orbital scenarios 

We consider two phenomenological scenarios based on the currently available 
MAGIC observations. The first one is such that there is an strict orbital periodicity 
of the neutrino emission, with the latter being correlated with the y-ray maximum 
measured by MAGIC (Albert et al. 2005), say, within a period of about 10 days 
(the active time) around phase 0.6 of the orbit. That is, the y-rays are produced 
during the time span where they have been detected and there is no production of 
y-rays at periastron. A scenario where this could happen even in a hadronic de- 
scription of the y-ray emission detected is briefly discussed below. A second case 
considered here is one in which there is orbital modulation of the y-ray emission, 
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Time [yrs] 



- Periodic signal 
Modulated signal 



Number of orbital periods 



Figure 1: Maximum neutrino to photon ratio compatible with the 95% CL interval for 
the corresponding neutrino background in the given integration time. See text for a 
more detailed explanation. 



i.e., y-rays are emitted along the orbit following the modulation imposed by the 
target matter and the absorbing photon fields. The maximum of the neutrino emis- 
sion would then be naturally expected at periastron, where the y-ray production 
should also be maximal but completely or nearly completely quenched by absorp- 
tion. This increases the active time along the orbit (we consider an active time of 
20 days) when significant neutrino emission proceeds. As a benchmark, we show 
in Table 1 the results of the neutrino yield for background and source events when 
there is no neutrino to photon enhancement (i.e., photons are not absorbed). In this 
situation, both for the cases of strict orbital periodicity and for the more relaxed 
modulated signal, the background greatly dominates the expected number of sig- 
nal events. The question we pose then is how much enhancement of the neutrino 
to photon ratio would still be compatible with a non-detection of the system in a 
km 3 -observatory. 



5.2.3 Numerical results 

To fix first the numerics given by the current constraints on neutrino emission, we 
consider that the neutrino flux follows the MAGIC spectrum and is enhanced up to 
the level constrained by AMANDA-II data (Ackerman 2006). In this situation, in 
the case of periodicity, the signal events for energies below 1 Te V in a detector such 
as ICECUBE during each orbital active time is 0.07 (with expected background, 
computed as explained above, of 0.12), whereas it is again 0.07 (but with back- 
ground 0.05) for energies above 1 TeV. These quantities maybe slightly increased 
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in a more comprehensive treatment of the detector, what we explore below. In the 
case of modulation, these quantities are (at least) a factor of 2 larger, given that 
the signal and background are proportional to the integration time. In one year 
we have 14 LS I +61 303 orbits, such that we expect about 2 signal events at all 
energies against 2.5 of background in the case of periodicity; or 4 signal events at 
all energies against a background of 5, in the case of modulation. Detecting, in 
just one year, only the expected number of background events or less, i.e., appar- 
ently no signal, in any of the cases mentioned, would then discard the theoretically 
predicted signal at the 68% CL (Feldman & Cousins 1998). 

Instead of fixing a priori the level of enhancement of the neutrino over the y-ray 
flux, it can be consider as a variable upon which constraints can be imposed from 
data. The 95% confidence level interval for the mean number of 1-year signal 
events in the case the number of detected events is equal to the expected back- 
ground is [0.00,5.75] and [0.00,6.26], for the periodicity and modulation cases, re- 
spectively. The amplification factor between the neutrino and photon fluxes is then 
limited to that which produces an expected signal that is smaller than the upper end 
of these intervals. We obtain ~18 and ~10, respectively. In one year, then, a neg- 
ative ICECUBE result for LS I +61 303 (i.e., a detection of a number of neutrino 
events that is compatible with the background expectation and are coming, within 
experimental uncertainty, from the direction of LS 1+61 303) would imply that the 
amplification factor between the neutrino and y-ray flux is severely constrained. In 
addition, the level of constraint will grow with additional observation time. Figure 
1 describes this improvement. As time goes by, the number of orbital periods that 
are integrated is larger. This increases the background and signal events expected 
from the theoretical assumptions. If just background events continue to be mea- 
sured, Feldman & Cousins' (1998) Table 4 and 5 gives the upper end of the 95% 
CL interval for the mean number of signal events that is still in agreement with the 
experimental data. Assuming that the neutrino flux is enhanced from that coming 
from the y-ray spectrum measured, the maximum neutrino to photon ratio that is 
compatible with the upper end of that 95% CL interval is obtained using Eq. J 1 3 b - 
This is what is plotted on the y-axis of the Figure 1 . After a few years of observa- 
tions of just the background neutrino events, the enhancement of the neutrino over 
the y-ray flux allowed is so low that theoretical models accounting for typical y-y 
absorption (e.g., see Christiansen et al. 2006) cannot accommodate it. Hadronic 
production of y-rays would thus be indirectly ruled out, with growing confidence. 
Alternatively, we can of course find the reverse result: the predicted (lower limit) 
signal from LS I +61 303 in muon neutrinos could lead to its discovery in ICE- 
CUBE. 

5.2.4 A more detailed treatment of the detector 

In Figure[2]we show the event rate above an energy threshold, shown on the hori- 
zontal axis using a more detailed numerical description of the ICECUBE detector. 
Specifically, in this (alternative) calculation we obtain the neutrino flux using the 
bolometric method (Alvarez-Muniz & Halzen 2002) that i) describes the spectra 
assuming the pionic origin of both gammas and neutrinos and ii) imposes equal 
energy in the flux of neutral and charged pions.We subsequently calculate the num- 
ber of Vf, induced neutrino tracks in IceCube using the semianalytical calculation 
presented by Gonzalez-Garcia, Halzen & Maltoni (2005) with quality cuts on the 
IceCube data referred to as level 2 cuts by Ahrens et al. (2003). This anticipated 



10 




Figure 2: Event rate of neutrinos above an energy threshold, with a full account of 
ICECUBE effective area. 



performance of the IceCube detector is consistent with the initial data from the first 
9 strings (see Achterberg et al. 2006). When the spectra have equal high energy 
slopes, this method will reproduce our previous results, i.e., can be seen that these 
results agrees well with the result already quoted above, where an effective area of 
a kilometer squared has been assumed. If we assume that the gamma-ray flux is 
steeper than the neutrino flux because of cascading, we get larger numbers for the 
neutrino predictions. These are also shown in the plot. 

5.2.5 Electron neutrinos 

A concomitant way of looking for neutrinos from LS I +61 303 would be the 
detection of electron neutrinos (see, e.g., Section IIB of Anchordoqui & Halzen 
2005 for details and formulae). The angular resolution of ICECUBE for the de- 
tection of such neutrinos is not yet known but expected to several degrees (say, 
several degrees) although with very large tails. That means that for a fraction of 
the events, say half, it may be similar to the AMANDA distribution: For those 
events with low uncertainty in angular position, the background increase because 
of the larger integration region would be compensated by the reduced atmospheric 
v e background (a factor of ~20 when compared with at 1 TeV), thus probably 
producing a similar confidence signal. The putative detection of LS I +61 303 in 
two different channels would provide further confirmation of the hadronic origin of 
the high energy radiation. Interesting to note in this case is that the acceptance in 
showers is 4n and ICECUBE would be able see the Southern micro-quasars as well 
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6 Discussion 



Note that for the former computations we have not made use of any particular 
model of the source or the nature of the system. Whatever they are, if the origin 
of the MAGIC-detected gamma-ray emission is hadronic, the neutrino to photon 
ratio can be constrained by ICECUBE, we have shown. We now briefly discuss our 
results in relation to published models of hadronic interpretations for LS I +61 303. 
Still, as we mentioned, a hadronic modelization of this system as a pulsar wind - 
stellar wind interaction is not available, so we focus on microquasar models. 

A model of y-ray and neutrino emission from microquasars based on hadronic 
interactions of relativistic protons in the jet with stellar wind ions have been pro- 
posed by Romero et al. (2003) and Torres et al. (2005). This kind of modelling has 
several free parameters and embedded assumptions, and in order to test its viabil- 
ity, it is of interest to compare the obtained results with the report of the MAGIC 
detection of LS I +61 303. The specific application of the quoted model to LS 
I +61 303 has been developed by Romero et al. (2005) and Christiansen et al. 
(2006). A prediction of this model as obtained in the previous works is that the 
maximum of the y-ray radiation is to be found at periastron (even when consider- 
ing in some detail the orbital dependence of y-y absorption, which is also maximal 
there). This is not supported by MAGIC results, which have only imposed upper 
limits at this orbital phase and found that the maximum of the emission is instead 
shifted to the phases where the radio and X-ray maxima are also located, after 
the periastron passage. Is this fact enough to rule out this particular hadronic in- 
terpretation (i.e., in the framework of the wind-jet interacting model) in favor of 
leptonic ones (e.g., Bednarek 2006)? In applying the hadronic modelling to LS I 
+61 303, and in order to motivate a hadronic interpretation, Romero et al. (2005) 
have emphasized the anti-correlation of the GeV and radio emission reported by 
Massi 2004. This, however, overestimates the goodness and orbital coverage of 
EGRET data on LS I +61 303. The typical EGRET viewing period has a similar 
duration than the timescale that is to be tested, i.e., the orbital period, which casts 
doubts upon the resulting anti-correlation even if all other caveats (e.g., smallness 
of the data sample and errors for each of the lightcurve points) were to be put aside. 
EGRET data can be used to claim variability of the emission in monthly timescales 
(Torres et al. 2001, Nolan et al. 2003), and even that is not conclusive due to the 
large error bars of each of the flux data points. Assigning a precise orbital modula- 
tion at the GeV regime is clearly a task for GLAST. The former satellite will have 
a sensitivity about 30-50 times better than EGRET, and in a matter of days will 
be able to detect the weakest of the EGRET sources. Recent simulations (Dubois 
2006) show that GLAST will be able to follow the variabilty of LS I +61 303 in 
intra-orbit integrations, demonstrating the ability to detect the period of the system, 
should the GeV emission be correlated (or anti-correlated) with it. Massi (2004) 
proposed that the EGRET maximum is correlated with the periastron of the sys- 
tem and Romero et al. and Christiansen et al. used this fact, requiring that the 
GeV flux at periastron is consistent with the EGRET measurement, in order to fix 
a free phenomenological parameter that accounts for some wind particles being 
unable to diffuse into the microquasar jet. This parameter (so-called /,,) linearly 
affects the computation of the y-ray (and neutrino) luminosity. f p was therefore 
chosen to be 0. 1 all along the orbit, this reduces one order of magnitude the other- 
wise predicted GeV and TeV luminosity. The parameter f p is ad-hoc, thus, it is 
not defined by any formula. It is likely that several effects can impede wind ions 
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to enter (diffuse) into the jet. Shock formation on the boundary layers or partial 
alignment between jet and wind, are examples. For instance, by comparing the 
diffusion and convection timescales it can be seen that efficient diffusion of ions 
into the jet is only possible when the wind blows from the side (Torres et al. 2004). 
Thus, if the wind-jet inclination changes with time, either by geometry or wind 
dynamics, the efficiency in the diffusion of particles will also change. But these 
effects may not be constant along the orbit. In fact, there is no reason to expect so 
when the wind velocity at the impact with the jet, the accretion rate, the distance 
between stars, and the wind density all significantly vary (order of magnitude) be- 
tween periastron (phase 0.23) and the y-ray maximum (phase ~ 0.6). In addition 
also the jet inclination with respect of the stellar wind flow, particularly if the latter 
is not strictly confined to the orbital plane, or if the angle of confinement change 
randomly, may change with time. Just phenomenologically, then, if the penetra- 
tion factor f p changes significantly from to 1 along the orbit, even in a hadronic 
model, both strict periodicity (as commented in the previous section), and thus, 
no radiation at periastron, and a modulated scenario with emission all along the 
orbit are possible. This is a conceptually new possibility for hadronic models in 
microquasars, which may even be useful for other cases should the nature of LS I 
+61 303 finally reveals as jetless, where radiation is quenched not because of the 
existence of large stellar fields and thus an increased opacity, but rather because of 
a modulation of the target matter. This reminds of collective stellar wind mod- 
els (Romero & Torres 2003, Domingo-Santamarfa & Torres 2006). The fact that 
the phases at which the y-ray maxima occurs is located at the second maximum 
of the accretion rate over the compact object provides the needed target matter 
for a hadronic interpretation of the y-ray detection at this orbital phase (provided 
the diffusion is efficient, f p ~ 1), whereas it reduces the resulting enhancement 
factor between the neutrino and y-ray flux (due to the reduced level of y-y opacity 
found far from the periastron of the system). This implies that longer integration 
times in neutrino telescopes such as ICECUBE are needed to reach ruling-out lev- 
els in the strict periodicity scenario (consistent with Fig. 1). Additional sources 
of opacity at periastron would decrease the level of y-ray flux, but increase the 
neutrino-to-photon ratio; ICECUBE will be able to test this possibility. Further 
MAGIC observations as well as ICECUBE results on LS 1+61 303, a comparison 
of these with the results of the previous section, and a more detailed implementa- 
tion of the conceptual ideas about suppression of target matter at periastron, will 
help further constrain, or finally rule out, a hadronic interpretation for this system. 
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3 Happening because the accretion rate is cc v~^, where v re i is the relative velocity between the compact 
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